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We have investigated macroscopic quantum tunneling (MQT) in in-plane high-T c superconductor 
Josephson junctions and the influence of the nodal-quasiparticle and the zero energy bound states 
(ZES) on MQT. We have shown that the presence of the ZES at the interface between the insulator 
and the superconductor leads to strong Ohmic quasiparticle dissipation. Therefore, the MQT rate is 
noticeably suppressed in comparison with the c-axis junctions in which ZES are completely absent. 

PACS numbers: 



INTRODUCTION 

A mesoscopic single Josephson junction is an inter- 
esting physical object for testing quantum mechanics at 
a macroscopic level. In current-biased Josephson junc- 
tions, the measurements of macroscopic quantum tun- 
neling (MQT) are performed by switching the junction 
from its metastable zero- voltage state to a non-zero volt- 
age state (see Fig. 1 (d)). Until now, experimental inves- 
tigations of MQT have been focused on s-wave (low-T c ) 
junctions only. This fact is due to the naive preconcep- 
tion that the existence of the low energy quasiparticles in 
the (i-wave order parameter of a high-T c cuprate super- 
conductor [l[ may preclude the possibility of observing 
the MQT. 

Recently we have theoretically investigated the effect 
of the nodal-quasiparticle on MQT in the d-wave c-axis 
junctions (e.g., Bi2212 intrinsic Josephson junctions 12|, 
hj and cross whisker junctions [lj|) 0, y]. We have 
shown that the effect of the nodal-quasiparticle gives rise 
to a super-Ohmic dissipation 0, [B| and the suppression 
of the MQT due to the nodal-quasiparticle is very weak. 

The first experimental observation of the MQT in the 
high-r c Josephson junction was made by Bauch et aL, us- 
ing a YBCO grain boundary bi-epitaxial junction [g, 0]. 
Recently, Inomata et al. @], Jin et al. and Kashiwaya 
et al. [IlJIlH have experimentally observed the MQT in 
the c-axis (Bi2212 intrinsic) junctions. They reported 
that the effect of the nodal-quasiparticle on the MQT 
is negligibly small and the thermal-to-quantum crossover 
temperature is relatively high (0.5~1K) compared with 
the case of the low-T c and the YBCO bi-epitaxial junc- 
tions. In Jin et alJ s experiment, 0(N 2 ) (N is the num- 
ber of the stacked junctions) enhancement of the MQT 
rate was reported. This enhancement is attributed to 



collective motion of the phase differences in the intrinsic 
junctions HI, 1(| 12 1. 

In this paper, we will theoretically investigate the 
MQT in the rf-wavc in-plane junctions parallel to the ab- 
plane (see Fig. 1) [ijj]. In such junctions, the zero energy 
bound states (ZES) [3] are formed near the interface 
between superconductor and the insulating barrier. ZES 
are generated by the combined effect of multiple Andreev 
reflections and the sign change of the (i-wave order pa- 
rameter symmetry, and are bound states for the quasi- 
particle at the Fermi energy. Below, we will show that 
ZES give rise to Ohmic type strong dissipation so MQT 
is considerably suppressed in compared with the c-axis 
and the do /do junction cases. 



(a) 



d 



^ d K 




W do 



U!) 



ZES ZES 




ZES 



FIG. 1: Schematics of the in-plane d-wave Josephson junc- 
tion, (a) do /do, (b) do/d„u, and (c) d^^/d^/4 junction. In 
the case of the do/d n /4 and d 7r /4/d n /4 junctions, the ZES are 
formed near the boundary between superconductor d^/4 and 
insulating barrier I. (d) Potential U((j>) v.s. the phase differ- 
ence 4> between two superconductors. Uo is the barrier height 
and Up is the Josephson plasma frequency. 
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EFFECTIVE ACTION 

By using the method developed by Eckern el. al., [2(| 
the partition function of the system can be described by 
a functional integral over the macroscopic variable (the 
phase difference </)), 



Z = 



J V^t) exp (- 



Seff[4>] 



(1) 



In the high barrier limit, i.e., z = mwo/h 2 kp 3> 1( m 
is the mass of the electron, wq is the height of the delta 
function potential I, and kp is the Fermi wave length), 
the effective action S e ff is given by 



S, 



e//L 



no 



dr 



Tui 



M fd(j){T) 
T \~dr~ 

hp 



U(</>) 



S a [$], 



S a [$>] = dr dr'a{T - r') cos 



(2) 

In this equation, j3 = 1/fcsT, M — C (h/2e) 2 is the mass 
(C is the capacitance of the junction) and the potential 
U(4>) is described by 
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where Ij is the Josephson current and I ex t is the external 
bias current, respectively. The dissipation kernel a(r) is 
related to the quasiparticle current I qp under constant 
bias voltage V by 
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(4) 



at zero temperature. 

Below, we will derive the effective action for the 
three types of the d-wave junction (do/do, do/d v u, and 
d-K/i/dn/4) in order to investigate the effect of the nodal- 
quasiparticles and ZES on MQT. In the case of the do /do 
junction, the node-to-node quasiparticle tunneling can 
contribute to the dissipative quasiparticle current. How- 
ever, ZES are completely absent. These behaviors are 
qualitatively identical with that for the c-axis Joseph- 
son junctions @, Cm the other hand, in the case of 
the do/d n /4 and d^/i/d^/^ junction, the ZES are formed 
around the surface of the superconductor d„u. There- 
fore the node to ZES (do/d n /i) and the ZES to ZES 
{d^/i/dir/i) quasiparticle tunneling becomes possible. 

Firstly, we will calculate the potential energy U in the 
effective action (2). As mentioned above, U can be de- 
scribed by the Josephson current through the junction in 
the high barrier limit. In order to obtain the Josephson 
current we start from the Bogoliubov-de Gennes (B-dG) 
equation [l9j . 
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where cp is the phase of order parameter, u(v) is the 
amplitude of the wave function for the electron (hole)- 
like quasiparticle, h(r) — — h 2 X7 2 /2m — \i + w Q S(x), and 

A(r — r') — SI -1 J2k ex P ^ ' ( r — r ')] ^ s * ne 0T< ^ eT 
parameter (f2 is the volume of the superconductor). In 
the superconductor regions (do and d^/4), the B-dG equa- 
tion (5) can be transformed into the eigenequation 



A 



k e 



A fe e* 



U k\ =E 
v k 



u k 
"k 
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where, £u = h 2 k 2 /2m + h 2 p 2 /2m — \i (p = 2nn/D and D 
is the width of the junction). The amplitude of the or- 
der parameter A^ is given by A cos 29 = A do (9) for d 
and Aosin20 = Ad w/i (9) for d^/4, where cos9 = k/kp. 
The Andreev reflection coefficient for the electron (hole)- 
like quasiparticle r^ e (r e /i) is calculated by solving the 
eigenequation (6) together with the appropriate bound- 
ary conditions. By substituting rh e {r e h) into the formula 
of the Josephson current for unconventional superconduc- 
tors (the Tanaka-Kashiwaya formula) [l9[ , 
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we can obtain <f> dependence of the Josephson current. 
Here, A± = A, 



n± = V(frtu n ) 2 - |A±| 



(2n + 1)tt/PH is the fermionic Matsubara frequency. In 



the case of low temperatures (/? 
barrier limit (zq 3> 1), we get 



<C Aq) and the high 
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where h = 3TrA /10eR N , h = TT 2 h(3Al/35e 3 N c R 2 N , and 
^3 = 37rzoAo/4ei?Ar (Rn = 3Trhz 2 /2e 2 N c is the normal 
state resistance of the junction and N c is the number of 
channel at the Fermi energy). 

By substituting the Josephson current into eq. (3), we 
can obtain the analytical expression of the potential U, 
i.e., 
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As in the case of the s-wave and the c-axis junctions 0], 
U can be expressed as the tilted washboard potential (see 
Fig. 1(d)). 
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DISSIPATION DUE TO 
NODAL-QUASIPARTICLES AND ZES 

Next we will calculate the dissipation kernel ot(r) in 
the effective action (2). In the hig h barrier limit, the 
quasiparticle current is given by [19( 



2e 



Em 2 f 

m J — ( 
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x [f(E) — f(E + eV)\ , 



(10) 



where tjy « cos9/zq is the transmission coefficient of the 
barrier I, N L ^(E,9) is the quasiparticle surface den- 
sity of states (L = do and R = do or d^/4) and f{E) is 
the Fermi-Dirac distribution function. The explicit ex- 
pression of the surface density of states was obtained by 
Matsumoto and Shiba [2l| . In the case of do , there are 



no ZES. Therefore the angle 9 dependence of Nd (E,t 
is the same as the bulk and is given by 
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On the other hand, iV c ^ 7r/4 (£ , , 9) is given by 
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The delta function peak at E = corresponds to the 
ZES. Because of the bound state at E = 0, the quasi- 
particle current for the do/d^/4 and d 7r /4/d 7r /4 junctions 
is drastically different from that for the do /do junctions 
in which no ZES are formed. By substituting Eqs. (11) 
and (12) into Eq. (10), we can obtain the analytical ex- 
pression of the quasiparticle current I qp (V). In the limit 
of low bias voltages (eV <C A ) and low temperatures 
(/3 _1 <C Aq), this can be approximated as 
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In the calculation of I qp for the d^/^/d^u junctions, we 
have replaced the ZES delta function 6(E) in Eq. (12) 
with the Lorentz type function, i.e., 
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in order to avoid a mathematical difficulty and model 
the real systems (which include e.g. disorder and many 
body effects). It is apparent from Eq. (13) that, in the 
case of do /do junctions, the dissipation is of the super- 
Ohmic type as in the case of the c-axis junctions Q. 



This can be attributed to the effect of the node-to-node 
quasiparticle tunneling. Thus the quasiparticle dissipa- 
tion is very weak. On the other hand, in the case of the 
do/dir/4 junctions, the node-to-ZES quasiparticle tunnel- 
ing gives the Ohmic dissipation which is similar to that 
in normal junctions [2(J. Therefore the dissipation for 
the do/d^/4 junctions is stronger than that for the do/ do 
junctions. Moreover, in the case of the d v 1 4 / d v / 4 junc- 
tions, the ZES-to-ZES quasiparticle tunneling dominates 
the quasiparticle dissipation. The broadening of the ZES 
peak e is typically one order magnitude smaller than Ao. 
Therefore, due to the prefactor (A /e) 2 in Eq. (12), the 
quasiparticle dissipation in the d n /4/d 7r /4 junctions be- 
comes enormously stronger than that for the do/ do and 
do/d-x/4 cases. 

From Eq. (4), the asymptotic form of the dissipation 
kernel is given by 
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The result for dp /dp junction is in agreement with previ- 
ous works di,E3,[2i|. 

MQT IN IN-PLANE d-WAVE JUNCTIONS 

Let us move to the calculation of the MQT rate T for 
the d-wave Josephson junctions based on the standard 
Caldeira and Leggett theory [24| . At zero temperature T 
is given by 



Aexpj-^ 



(16) 



where Sb = S e ff[4>B\ and 4>b is the bounce solution. 
Following the procedures as above, we obtain the ana- 
lytical formulae of the MQT rate for the in-plane d-wave 
junctions as 
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where c = f Q °°dy y 4 \n(l + 1/y 2 )/ sinh 2 (7n/) w 0.0135, 
£(3) is the Riemann zeta function, r\ = Rq/ R^ is the dis- 
sipation parameter, Up is the barrier height of the poten- 
tial U, uj p is the Josephson plasma frequency, x = I e xt/h 
(i=l,2,3), and 
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is the MQT rate without the dissipation. In Eq. (17) 



6AI 



24^2 A y_! 



dy y 



VI - y Jo 



dz z 2 K x (z\y\Y 



(19) 



is the renormalized mass due to the high frequency com- 
ponents (lo > Up) of the quasiparticle dissipation. 

In order to compare the influence of the ZES and the 
nodal-quasiparticle on the MQT more clearly, we will 
estimate the MQT rate (17) numerically. For a meso- 
scopic bicrystal YBCO Josephson junction [25| (Ao = 
17.8 meV, C = 20 x lO" 15 F, R N = 100 Q, x = 0.95), 
the MQT rate is estimated as 
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As expected, the node-to-ZES and ZES-to-ZES quasi- 
particle tunneling in the do/d^u and d^/n/d^/i junc- 
tions give strong suppression of the MQT rate compared 
with the do / do junction cases. Moreover in the d^/i/d^/i 
cases, the MQT is almost completely depressed. 



SUMMARY 

In conclusion, MQT in in-plane high-T c superconduc- 
tors has been theoretically investigate and analytically 
obtained the formulae of the MQT rate which can be 
used to analyze experiments. The node-to-node quasi- 
particle tunneling in the do /do junctions gives rise to the 
weak super-Ohmic dissipation as in the case of the c-axis 
junctions [2j. For the do/d^/4 junctions, on the other 
hand, we have found that the node-to-ZES quasiparti- 
cle tunneling leads to the Ohmic dissipation. Moreover, 
in the case of the d^/i/d^u junctions, the ZES-to-ZES 
quasiparticle tunneling gives very strong Ohmic dissipa- 
tion so the MQT is drastically suppressed. 

In this paper we have considered the high barrier limit 
case (zq 3> 1) only. In the low barrier cases, the ZES 
becomes split into two finite energy Andreev levels due 
to the ZES resonance Moreover, the energy of the 
split Andreev levels depends on the phase difference <j> 
and the influence of the proximity effect becomes more 
important. To take into account such effects, the present 
approach should be considerably modified. This issue 
will be investigated in future articles. 
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